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oo ; Abstract 

Dispersion relations provide a powerful tool to analyse the electromagnetic form factors of the 
nucleon for all momentum transfers. Constraints from meson- nucleon scattering data, unitarity, 
and perturbative QCD can be included in a straightforward way. In particular, we include the 
Ph! 2tt, pir, and KK continua as independent input in our analysis and provide an error band for our 

results. Moreover, we discuss two different methods to include the asymptotic constraints from 
perturbative QCD. We simultaneously analyze the world data for all four form factors in both 
the space-like and time-like regions and generally find good agreement with the data. We also 
extract the nucleon radii and the ujNN coupling constants. For the radii, we generally find good 
agreement with other determinations with the exception of the electric charge radius of the proton 
which comes out smaller. The uNN vector coupling constant is determined relatively well by the 
fits, but for the tensor coupling constant even the sign can not be determined. 
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I. INTRODUCTION 



The electromagnetic (em) form factors of the nucleon provide an important tool to study 
strong interaction dynamics over a wide range of momentum transfers 0, . Their detailed 
understanding is important for unraveling aspects of perturbative and nonperturbative nu- 
cleon structure. At small momentum transfers, they are determined by gross properties of 
the nucleon like the charge and magnetic moment. At high momentum transfers they encode 
information on the quark substructure of the nucleon as described by perturbative QCD. The 
form factors also contain important information on nucleon radii and vector meson coupling 
constants. Moreover, they are an important ingredient in a wide range of experiments from 
Lamb shift measurements to determinations of the strangeness content of the nucleon 0| . 

With the advent of the new continuous beam electron accelerators such as CEBAF (Jeffer- 
son Laboratory), ELSA (Bonn), and MAMI (Mainz), a wealth of precise data for space-like 
momentum transfers has become available p. Due to the difficulty of the experiments, the 
time-like form factors are less well known. While there is a fair amount of information on the 
proton time-like form factors ja l7[ , only one measurement of the neutron form factor from 
the pioneering FENICE experiment exists. Recently, new precise data on the proton 

time-like form factors have been presented by the BES, CLEO, and BABAR collaborations 

IQEil. 



In this work, we therefore analyze the nucleon form factors in both the space- and time- 
like regions. We use dispersion relations which provide a model-independent framework to 
consistently analyze the form factor data in both regions. A complete description of our 
data base is given in section |V| 

It has been known for a long time that the pion plays an important role in the long- 
range structure of the nucleon [l3j. This connection was made more precise using dispersion 
theory in the 1950's 0,0]. Subsequently, Frazer and Fulco have written down partial wave 
dispersion relations that relate the nucleon electromagnetic structure to pion-nucleon (ttN) 
scattering and predicted the existence of the p resonance 0, 0] • Hohler et al. (Hf first 
performed a dispersion analysis of the electromagnetic form factors of the nucleon including 
the 2n continuum derived from the pion form factor and 7rA^-scattering data llfll . In the 
mid-nineties, this analysis has been updated by Mergell, Meifiner, and Drechsel |20| and was 
later extended to include data in the time-like region MS- The new precise data for the 
neutron electric form factor have been included as well |23| . For recent attempts to calculate 
the nucleon form factors in QCD using nonperturbative methods, see e.g. Refs. jiil EiL 126^. 

A recent form factor analysis by Friedrich and Walcher [2jj created some renewed interest 
in the 2n continuum. They analysed the electromagnetic nucleon form factors and performed 
various phenomenological fits [27| . Their fits showed a pronounced bump-dip structure in 
G\ which they interpreted as a signature of a very long-range contribution of the pion cloud 
to the charge distribution in the Breit frame extending out to about 2 fm. This observation 
is at variance with the pion cloud contribution to the nucleon form factors as given by the 27r 
continuum - the lowest-mass intermediate state including pions only 0| . We will discuss a 
possibility how to reconcile these findings below. 

It is well known that vector meso ns p lay an important role in the electromagnetic struc- 
ture of the nucleon, see e.g. Refs. 000 El ESQ, and the remaining; contributions 



to the spectral function have usually been approximated by vector meson resonances. A novel 



i py 

addition in this work is the inclusion of the KK |35l . l36j | and pir |37( continuum contribu- 
tions similar to the 2tt continuum described above. The continuum contributions provide 
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independent information on the spectral functions that supplements the electromagnetic 
form factor data. Moreover, we enforce the asymptotic constraints from pQCD. We discuss 
two different approaches to obtain this behavior: super convergence relations and an explicit 
continuum term with the correct pQCD behavior, and show results for both methods. 

Our manuscript is organized as follows. Sec. Ill Al contains the basic definitions of the 
nucleon em form factors and a short discussion of the corresponding dispersion relations. 
The various continuum contributions to the spectral functions are discussed in Sec. 1111 Al 
The structure of and the constraints on the spectral functions are given in Sec. IIVI The 
results of our fits are presented and discussed in Sec El We close with a short summary and 
outlook in Sec. IVII The fit parameters are collected in the appendix. 



II. PRELIMINARIES 



A. Definitions 



The electromagnetic (em) structure of the nucleon is determined by the matrix element 
of the vector current operator between nucleon states as illustrated in Fig. [T] 




FIG. 1: The nucleon matrix element of the electromagnetic current j^ m . 

Using Lorentz and gauge invariance, this matrix element can be expressed in terms of 
two form factors, 



{P'\37\p)=u{p') 



. , F 2 (t) 



u(p) 



(1) 



where M is the nucleon mass and t = (p'—p) 2 the four- momentum transfer squared. For data 
in the space-like region, it is often convenient to use the variable Q 2 = —t > 0. The functions 
Fi(t) and F 2 (t) are the Dirac and Pauli form factors, respectively. They are normalized at 
t = as 

Ff (0) = 1 , FT (0) = , Ff (0) = k p , F,T(0) = K n , (2) 

with K p = 1.793 and K n = —1.913 the anomalous magnetic moment of the proton and the 
neutron, respectively, in units of nuclear magnetons. 

It is convenient to work in the isospin basis and to decompose the form factors into 
isoscalar and isovector parts, 



(3) 



where i = 1,2. The experimental data are usually given for the Sachs form factors 

G E (t) = F 1 (t)-TF 2 (t), 

Gm if) = F 1 (t) + F 2 (t), (4) 
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where r = — t/(4M 2 ). In the Breit frame, Ge and Gm may be interpreted as the Fourier 
transforms of the charge and magnetization distributions, respectively. 

The nucleon radii a/ (r 2 ) can be defined from the \ow-t expansion of the form factors, 

F(t) = F(0)[l+t(r 2 )/6 + ...] , (5) 

where F(t) is a generic form factor. In the case of the electric and Dirac form factors of the 
neutron, G E and F™, the expansion starts with the term linear in t and the normalization 
factor F(0) is dropped. 



B. Dispersion Relations and Spectral Decomposition 

Based on unitarity and analyticity, dispersion relations relate the real and imaginary 
parts of the electromagnetic nucleon form factors. Let F(t) be a generic symbol for any one 
of the four independent nucleon form factors. We write down an unsubtracted dispersion 
relation of the form 

1 f 30 ImF(t') , , 
F(t) = - / - / { dt', 6 
k J t t'-t-ie 

where to is the threshold of the lowest cut of F(t) (see below) and the it defines the integral 
for values of t on the cut. The convergence of an unsubtracted dispersion relation for the 
form factors has been assumed. We could also use a once subtracted dispersion relation, since 
the normalization of the form factors at t = is known. Using Eq. © the electromagnetic 
structure of the nucleon can be related to its absorptive behavior. 

The imaginary part Im F entering Eq. (JBj) can be obtained from a spectral decomposition 
0, For this purpose it is most convenient to consider the electromagnetic current 

matrix element in the time-like region (t > 0), which is related to the space-like region 
(t < 0) via crossing symmetry. The matrix element can be expressed as 

Jfi = (N(p)N(p)\jT(m (7) 



u(p) 



v(p) 



cm 



where p and p are the momenta of the nucleon and antinucleon created by the current j . 
respectively. The four-momentum transfer squared in the time-like region is t — (p + p) 2 . 

Using the LSZ reduction formalism, the imaginar y p art of the form factors is obtained 
by inserting a complete set of intermediate states as jlj, HH 

ImJ, = |(27r) 3 /W^(p|J w (0)|A)(A|j-(0)|0)^(p)5 4 (p + p-p A ), (8) 

A 

where Af is a nucleon spinor normalization factor, Z is the nucleon wave function renormal- 
ization, and Jn{%) = J ( x )lo with Jn(x) a nucleon source. This decomposition is illustrated 
in Fig. EJ It relates the spectral function to on-shell matrix elements of other processes. 

The states |A) are asymptotic states of momentum p\. They carry the same quantum 
numbers as the current j^ m : I G (J PC ) = 0~(1 ) for the isoscalar current and I G (J PC ) = 
1 + (1 ) for the isovector component of j^ m . Furthermore, they have zero net baryon number. 
Because of G-parity, states with an odd number of pions only contribute to the isoscalar 
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FIG. 2: The spectral decomposition of the nucleoli matrix element of the electromagnetic current 
j e ™. | A) denotes an hadronic intermediate state. 



part, while states with an even number contribute to the isovector part. For the isoscalar 
part the lowest mass states are: 3tt, 5tt, . . ., KK, KKtc, . . .; for the isovector part they are: 

27T, 47T, . . .. 

Associated with each intermediate state is a cut starting at the corresponding threshold 
in t and running to infinity. As a consequence, the spectral function lm.F(t) is different 
from zero along the cut from t to oo, with i = 4 (9) M% for the isovector (isoscalar) case. 

The spectral functions are the central quantities in the dispersion-theoretical approach. 
Using Eqs. (J7I8)) . they can in principle be constructed from experimental data. In practice, 
this program can only be carried out for the lightest two-particle intermediate states. 

The longest-range, and therefore at low momentum transfer most important continuum 
contribution comes from the 2ti intermediate state which contributes to the isovector form 
factors . A new calculation of this contribution has recently been performed in Ref . IS . 
In this analysis we for the first time also include the KK and pir continua 



aa Balsa- 



III. CONTINUUM CONTRIBUTIONS 

Our general strategy is to include as much physics information in the construction of 
the spectral functions as possible. In this section, we explicitly contruct the 27r, KK, and 
pit continua mentioned above. These continua will be an important part of our spectral 
functions. 



A. 2ir Continuum 

The 27r contribution has recently been reevaluated in a model-independent way |38[ using 
the latest experimental data for the pion form factor from CMD-2 @, KLOE 0], and SND 



41[. Here we give a short summary of this evaluation. 

Following Refs. H, 4^|, the 2n contribution to the isovector spectral functions in terms 



of the pion charge form factor F n (t) and the P-wave nn — > NN amplitudes f±(t) can be 

1 Note that the effect of these continua has previously been studied in Ref. |3^| using fits to parameterizations 
of form factor data. 
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expressed as: 



lmG° u (t) = $F„{t)*f l _{t) , (9) 



where q t = y/t/A — M 2 . The imaginary parts of the Dirac and Pauli Form factors can be 
obtained using Eq. (JIJ). The 2n continuum is expected to be the dominant contribution to 
the isovector spectral function from threshold up to masses of about \ft ~ 1 GeV 4^ • 



The P-wave -kit — > NN amplitudes flit) are tabulated in Ref. |42j]. The representation 
of Eq. Q gives the exact isovector spectral functions for 4M 2 <t< 16M 2 , but in practice 
holds up to t ~ 50M 2 « 1 GeV 2 . Since the contributions from 4tt and higher intermediate 
states is small up to t ~ 50M 2 , F n (t) and the f±(t) share the same phase in this region and 
the two quantities can be replaced by their absolute values. 

The experimental data for the pion form factor from CMD-2 Q, KLOE H}, and SND 



41 



show some discrepancies. In Ref. [38], the 2tt continuum given by Eq. © was evaluated 
for all three sets and the errors from the discrepancy between the sets were estimated. The 
resulting difference in the spectral functions is very small 1%). It is largest in the p-peak 
region, but this region is suppressed by the ttti — ► NN amplitudes f±(t) which show a strong 
fall-off as t increases. 

The spectral functions have two distinct features. First, as already pointed out in (l6| . 
they contain the important contribution of the p-meson with its peak at t ~ 30M 2 . Second, 
on the left shoulder of the p, the isovector spectral functions display a very pronounced 
enhancement close to the two-pion threshold. This is due to the logarithmic singularity 
on the second Riemann sheet located at t c = 4M 2 — M^/M 2 = 3.98M 2 , very close to the 
threshold. If one were to neglect this important unitarity correction, one would severely 
underestimate the nucleon isovector radii {43j. In fact, precisely the same effect is obtained 
at leading one-loop accuracy in relativistic chiral perturbation theory [.This topic 

was also discussed in heavy baryon chiral perturbation theory (ChPT) |46l . |47j and in a 
covariant calculation based on infrared regularization |48l . |49|. Thus, the most important 
2tt contribution to the nucleon form factors can be determined by using either unitarity or 
ChPT (in the latter case, of course, the p contribution is not included). 

The contribution to the nucleon form factors is obtained by inserting the 2tt contribution 
to the spectral function into the dispersion relations Eq. (|B|). The result can be parameterized 
as 

F (vM m _ ai + bi(l-t/a)- 2 / i 

^ {t) ~ — 2(1 — ' '- 1 ' 2 ' (10) 

where a x = 1.10788, &i = 0.109364, c x = 0.36963 GeV 2 , d x = 0.553034 GeV 2 , a 2 = 5.724253, 
b 2 = 1.111128, c 2 = 0.27175 GeV 2 , and d 2 = 0.611258 GeV 2 . The errors in these constants 
are of the order 4% or less. Of course, it is not necessary to use a parameterization like 
Eq. (|10|) but the numerical evaluation of the fits becomes much simpler. The form factor 
contributions from Eq. (J10|) are shown in Fig. |21 below. 
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B. KK Continuum 



The KK contribution to the isoscalar spectral function was evaluated in Refs. [Til . 35l . 3 
from an analytic continuation of KN scattering data. In the following, we give a short 
summary of this work. 

The KK contribution to the imaginary part of the isocalar form factors is given by 3^, 36| 



luiF[ s ' KR \t) = Re 



lmF ( 2 s ' KR \t) = Re 



MqA 


' y/t 


*fi) 


.2V2M 


Mq t \ 


6 iAi/2 







b\ /2 '- 1/2 (t)-b{ /2 ' 1/2 (t) 



F K {ty , (11) 



F K {ty , (12) 



with p t = \JtjA — M 2 and qt = yt/A — Mg. Fx(t) represents the kaon form factor whereas 

the b\^ 2,±1 ^ 2 (t) are the J = 1 partial wave amplitudes for KK — > NN [35|,|36j]. Once these 
imaginary parts are determined, the contribution of the i^.K'-continuum to the form factors 
is obtained from the dispersion relation Eq. (JBJ). 

For t > AM 2 the partial waves are bounded by unitarity, 

Vpjqt\b\ /2 ' ±1/2 (t)\<l. (13) 

In the unphysical region AM\ <t< AM 2 , however, they are not constrained by unitarity. In 

Ref. [3.5], the amplitudes 6^ //2 ' ±1//2 (t) in the unphysical region have been determined from an 
analytic continuation of K iV-scattering amplitudes. The contribution of the physical region 
t > AM 2 in the dispersion integral © is suppressed for small momentum transfers and 
bounded because of Eq. |T3|) . Using the analytically continued amplitudes in the unphysical 
region and the unitarity bound in the physical region, the contribution of the KK continuum 
can therefore be calculated. Strictly speaking this calculation provides an upper bound on 
the spectral function since we replace the amplitudes and the form factor in Eqs. (fTT| IT2*|) 
by their absolute values. 

The striking feature in the spectral function is a clear <fr resonance structure just above 

— 1/2 1/2 

the KK threshold. The resonance structure appears in the partial wave amplitude b± ' 
as well as in the kaon form factor Fx- In contrast to the 2ir continuum, there is no strong 
enhancement on the left wing of the <ft resonance which sits directly at the KK threshold. 

The resulting contribution to the nucleon form factors can be parameterized by a pole 
term at the <fi mass: 

,{s,kk) ( ^_L I im^ (*J^/_ 0-i 



F ^* )(t) = _ yL± dt >~_£ 2 = 1,2, (14) 

V ; vr J AM 2 t'-t Ml-V 



K 



with dy K = 0.1054 GeV 2 and a^ K = 0.2284 GeV 2 . As a consequence, the contribution of 
the KK continuum to the electromagnetic nucleon form factors can conveniently be included 
in the analysis via Eq. (JHJ). The form factor contributions from Eq. (fT4^) are also shown in 
Fig. El below. 



C. pir Continuum 

Drawing upon a realistic treatment of the correlated pir exchange in the Bonn-Jiilich NN 
model , the pn contribution to the isoscalar spectral function was calculated in Ref. j3?J . 
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The contribution of the pn continuum can be evaluated in terms of a dispersion integral 
which in turn can be represented by an effective pole term for a fictitious u' meson with a 
mass Mv = 1.12 GeV El: 



F 



(t) 



ImK 



t'-t 



dt' 



M 2 ,,-t 



1,2 



(15) 



with af*' = —1.01 GeV 2 and af^ = —0.04 GeV 2 . In our form factor analysis, we use this 
effective pole instead of the full spectral function. 

There is very little sensitivity in our fits to , which can vary between —0.04 and —0.4 
without affecting the outcome of the fit. If the uj' pole is treated as a real resonance, the 
latter value is consistent with f u i ~ 10 for = —1.01 if the coupling constants g l w i NN 
(i = 1,2) from Ref. [33] are used as input. 

In Fig. El we show the contribution of the 27r, KK, and pir continua to the electromagnetic 




Q 2 [GeV 2 ] 



Q 2 [GeV 2 ] 



FIG. 3: The continuum contributions to the nucleon form factors F\ (left panel) and F2 (right 
panel) in the space-like region. The contribution of the 2tt continuum to the isovector form factors 
is given by the solid line, while the contribution of the KK and pir continua to the isoscalar form 
factors are given by the dash-dotted and dashed lines, respectively. 



nucleon form factors F% and F2. The 2ir contributes to the isovector form factors while the 
KK and pre continua contribute to the isoscalar form factors. The KK and pir contributions 
have opposite sign and partially cancel each other. The dominant contribution to F* comes 
from the pir continuum while for F| the KK contribution is larger. While the KK and pn 
contributions can be represented by simple pole terms, the expressions for the 2-7T continuum 
Eq. ()10j) are somewhat more complicated. This is related to the strong enhancement close 
to the 27r threshold on the left wing of the p resonance discussed above. Finally, note that 
these continuum contributions enter as an independent input in our analysis. They are not 
fitted to the form factor data. 
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IV. SPECTRAL FUNCTIONS 



A. Structure 

As discussed above, the spectral function can at present only be obtained from unitarity 
arguments and experimental data for the lightest two-particle intermediate states (27r and 
KK) 0, 135L The pn continuum contribution has been calculated in the Bonn-Jiilich 
NN modeljU^ 

The remaining contributions to the spectral function can be parameterized by vector 
meson poles. On one hand, the lower mass poles can be identified with physical vector mesons 
such as the uj and the </>. In the the case of the 3tt continuum, e.g., it has been shown in ChPT 
that the nonresonant contribution is very small and the spectral function is dominated by the 
uj ji^l . The higher mass poles on the other hand, are simply an effective way to parameterize 
higher mass strength in the spectral function. Different parameterizations are possible and 
an explicit example will be discussed below in relation to the pQCD behavior. 

In all our fits the spectral function includes the 2ir, KK, and pn continua from unitarity 
and the uj pole. Note that we also include a pole at the <fi mass to account for explicit 
strength not included in the KK and pn continua. In addition to that there are a number 
of effective poles at higher momentum transfers in the isoscalar (si,S2, ...) and isovector 
channels (t>i,t>2, ...). The spectral function has the general structure 

1m F t '(t) = lmF^ KK \t)+lmF^ /m \t)+ ^ na{5{M$-t), z = l,2, (16) 

V=U),4>,S\,... 

Jm F?(t) = lmF^\t) + £ iraY6(M$ - t) , i = 1, 2 . (17) 

V=V!,... 

The masses of the effective poles are fitted to the form factor data. We generally do not 
include widths for the effective poles. However, in some of the fits we allow a large width for 
the highest mass effective pole in order to mimick the imaginary part of the form factors in 
the time-like region. We have performed various fits with different numbers of effective poles 
and including/excluding some of the continuum contributions. In Sec. EJ we will discuss 
the results of these efforts. 



B. Constraints 



The number of parameters in the fit function is reduced by enforcing various constraints. 
The first set of constraints concerns the low-t behavior of the form factors: We enforce the 
correct normalization of the form factors as given in Eq. (J2J). The nucleon radii, however, 
are not included as a constraint. In some earlier fits, we had also constrained the neutron 
charge radius to the value from low-energy neutron-atom scattering experiments 51, 5^. In 
the fits discussed below, this constraint is dropped since the fit value is compatible with the 
empirical range from Refs. |5jJ, |52j . 

Perturbative QCD (pQCD) constrains the behavior of the nucleon electromagnetic form 



factors for large momentum transfer. Brodsky and Lepage [53j find for Q — > 00 



Fi(t) 



Q 2 ( 



i+l) 



, l Q2 

]n[ Ql 



1,2, 



(18) 
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where Qn — Aqcd- The anomalous dimension 7^2 depends weakly on the number of 
flavors |53[. The asymptotic behavior of the form factors has recently also been studied 
in connection to the unexpected behavior of the ratio Q 2 : F 2 (Q 2 ) / ' Fi{Q 2 ) for the proton 
measured at Jefferson Lab and different expressions for the logarithmic corrections were 
found HHj]. (For a further discussion of the asymptotic behavior of the nucleon form 
factors for large space- like and time-like momenta, see Ref. 56[.) In the current analysis we 



implement only the leading power behavior of the form factors and details of the logarithmic 
correction are not relevant. Note that the logarithmic term in Eq. (j!8|) was included in some 
of our earlier analyses 2(1 21, |23| but had little impact on the fit. The particular way this 



constraint was implemented, however, lead to an unphysical logarithmic singularity of the 
form factors in the time-like region which we want to avoid in the current analysis. 

The power behavior of the form factors at large t can be easily understood from perturba- 
tive gluon exchange. In order to distribute the momentum transfer from the virtual photon 
to all three quarks in the nucleon, at least two massless gluons have to be exchanged. Since 
each of the gluons has a propagator ~ 1/t, the form factor has to fall off as 1/t 2 . In the 
case of F 2 , there is additional suppression by 1/t since a quark spin has to be flipped. The 
power behavior of the form factors leads to super convergence relations of the form 



lmFi(t)t n dt = 0, t = 1,2, (19) 

to 

with n = for F\ and n = 0, 1 for F 2 . 

The pQCD power behavior can be enforced in various ways. To obtain some information 
about the induced theoretical uncertainty, we discuss two different methods in more detail: 

1. Super convergence (SC) approach: 

The asymptotic behavior of Eq. (JTSj) is obtained by choosing the residues of the vector 
meson pole terms such that the leading terms in the 1 / t-expansion cancel. This leads 
to a spectral function consistent with the superconvergence relations Eq. (|T9*|) and the 
asymptotic behavior Eq. 1)180 . This method is similar to what was used in earlier 
works [20|, |2l|, |23| . Here we add a very broad resonance of the structure 



-,(1, broad),,, _ a i (Mf - t) 

{M 2 -t) 2 + T 2 



Fr oa ^(t) = ^ l:^ , <=i,2, i= s ,v, (20) 



in both the isovector and isoscalar form factors. The resonance parameterises con- 
tinuum contributions in addition to the 2tt, KK, and pit continua and generates an 
imaginary part of the form factors in the time-like region for t > 4M 2 . The residue 
a\ as well as the mass and width parameters Mi and Tj are fit to the data. The 
width parameter Tj is of the same order of magnitude as the mass Mj and comes out 
typically of the order of a few GeV in our fits. 

2. Explicit pQCD continuum approach: 

In addition to satisfying the superconvergence relations, Eq. ()19|). a term of the form 

F (I, P QCD) _ of i - 1 2 I-s V (21) 

*i - 1 _ c 2 t + b 2(_ t y+i, 1-S,V, [2L) 

which explicitly enforces the pQCD behavior, Eq. (fTHj) . is added to the fit function. 
Such a term behaves like an effective resonance pole for small values of t, and restores 
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pQCD behaviour explicitly at high values of t. The super convergence relations cancel 
the leading order terms in the 1/t expansion. This explicit pQCD term is consistent 
with a nonvanishing imaginary part of the form factors in the time-like region. Note 
that the parameters \ and c* are the same in the isoscalar and isovector channels 
while the residue a\ depends on the channel. This method allows for a smoother 
interpolation between the low-t and large-t regions compared to the SC approach. 
Because of this feature, one might expect obtaining fits with fewer parameters. 

The number of effective poles in Eqs. (1161 [T7j) is determined by the stability criterion 
discussed in detail in j5?J. In short, we take the minimum number of poles necessary to 
fit the data. The number of free parameters is strongly reduced by the various constraints 
(unitarity, normalizations, super convergence relations). More details will be given together 
with the fits in the next section. 



V. FIT RESULTS 

The fits have been performed using the Fletcher- Reeves and the Polak-Ribiere conjugate 
gradient algorithms implemented in the GSL library 58]. To ensure initial convergence sta- 
bility, Monte-Carlo sampling was performed over the whole physically acceptable parametric 
volume to obtain a number of parameter sets with acceptable starting x 2 values. 

The constraints dictated by the normalization and the pQCD conditions have been rep- 
resented in terms of a set of linear equations for the resonance residua. The equations are 
solved each iteration using the LU decomposition. 

Soft constraints on composite variables which depend on a set of fit parameters allow 
to impose an exponential well for the set of parameters as a whole, limiting deviations of 
the composite variable from its desired central value. These constraints are implemented as 
additive \ 2 terms of the general form 

^ = p [ x _ ( x )]2 _ (^2) ( 22) 

where (x) is the desired central value, and p is the constraint strength parameter which 
allows to stabilize fit convergence over the whole range of iterations adaptively, regulating 
the steepness of the exp onential well. This method was used in some earlier fits to constrain 
the neutron radius [59(. 

The error bands are obtained by allowing the total x 2 /dof of the fit to be in the interval 
[xiiiniXmin + ^X 2 } where Xmin * s the X 2 value of the best fit, and 5\ 2 is obtained from the 
lcr confidence interval p-value equations, 5\ 2 — 1-04. 

The data basis used in the fits is taken from Ref. j2?| and in addition also includes the 
new data that have appeared since 2003 and the time-like data 0, H, H, EH EH Ell EH EH 
M, El EI EH 1 EHlsil, | 6^, 1 7(j . The CLAS collaboration at Jefferson Lab has performed 
measurements of the neutron magnetic form factor for momentum transfers 0.6 < Q 2 < 
5 GeV 2 jzH]. These data are still preliminary and are therefore generally not included in our 



fits. In subsection IV Dl however, we discuss a fit where these preliminary data are included. 
The results for G^, G P E , G P M are normalized to the phenomenological dipole fit: 

G D (Q 2 )=(l + ^ f ) 2 , (23) 



m 



D 
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where m 2 D = 0.71 GeV 2 . 

Additionally, certain features of the form factor behaviour in specific Q 2 ranges can be 
enhanced during the fitting procedure by artificially decreasing the errors on the experimen- 
tal data in that region as seen by the fit. This allows to, for example, explore the conditions 
necessary to produce a pronounced bump-dip structure in G%, discussed in Sec. IV El 

In the time-like region, the neutron data do not participate in the fit as they are obtained 
from a single experiment. They are therefore a genuine prediction. 

A. Superconvergence (SC) Approach 




FIG. 4: The nucleon electromagnetic form factors for space- like momentum transfer in the SC 
approach. The solid line gives our best fit while the dashed lines indicate the error band obtained 
by the la deviation as discussed above. 

In Fig.HJ we show the results in the SC approach for all four form factors in the space-like 
region compared to the world data. In general, we get a good description of all data within 
our error bands. For G P E there is some inconsistency in the data points around Q 2 ~ 1 GeV 2 . 
Our fit favors the lower data points in this region. As can be seen from the inset, the data 
at low momentum transfers are well described. In addition to the u and the residual <f), this 
fit has 2 more isoscalar poles (M S1 rs 1.1 GeV, M S2 »s 2.0 GeV) and 5 isovector poles with 
masses ranging from 1 to 3 GeV. The heaviest poles in both channels are broad resonances 
(cf. Eq. (120(0 with width parameters ranging from 5 GeV (isoscalar) to 19 GeV (isovector). 
All other poles have zero widths. The fit has 17 free parameters and a total x 2 /dof of 1.8. 
The fit parameters are listed in detail in Table HP in the Appendix. 
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Note also that we do not obtain a pronounced bump-dip structure in G% as observed in 
Ref. J2?J. However, all data for G 1 ^ are described within our error band and the experimental 
error. We will come back to this bump-dip structure in Sec. IV El and discuss the modifications 
in the spectral function required to produce this structure. 

In Fig. EJ we show our fit results for the SC approach in the time-like region. As in the 




FIG. 5: The nucleon electromagnetic form factors for time-like momentum transfer in the SC 
approach (right panel: G^f , left panels: G^ ). The proton data participate in the fit while the 
neutron data are a genuine prediction. The solid line gives our best fit while the dashed lines 
indicate the error band obtained by the lcr deviation as discussed above. 

space-like region we get a good description of the world data within our error bands. Our 
best fit, however, cannot reproduce the strong rise of the data for G P M near threshold. The 
data for G 1 ^ are also well described. Note that the neutron data do not participate in the 
fit and the corresponding curves are therefore genuine prediction of the dispersion analysis 
based on data in the other channels. 



B. Explicit pQCD Continuum Approach 

In Fig. El we show our results in the explicit pQCD continuum approach for space-like 
momentum transfers compared to the world data. Again, we get a good description of all 
data within our error bands. In contrast to the superconvergence approach, our fit now favors 
somewhat larger values of G V E in the region around Q 2 ~ 1 GeV 2 . For G\ the situation is the 
same as before: we describe all data within our error band and the experimental errors but 
see no pronounced bump-dip structure in the fits. In addition to the u and the residual 0, 
this fit has one more isoscalar pole (M 81 ~ 1.8 GeV) and 3 isovector poles (M V1 ph 1.0 GeV, 
M V2 ph 1.6 GeV, and M V3 pa 1.8 GeV). Moreover, it contains an explicit pQCD continuum 
term, Eq. ([21)1 . as discussed above. This fit has 14 free parameters and a total x 2 /dof of 2.0. 
In general, we cannot get a satisfactory description of all data using fewer parameters. The 
parameters of this fit are given in detail in Tables UTTl and HVl in the Appendix. 

In Fig. [7[ we show the results in the time-like region. As before we describe the time-like 
data within our error band, but the errors increase strongly close to threshold. Our best fit 
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FIG. 6: The nucleon electromagnetic form factors for space- like momentum transfer with the 
explicit pQCD continuum. The solid line gives our best fit while the dashed lines indicate the error 
band obtained by the la deviation as discussed above. 



turns over very close to threshold and cannot describe the two lowest data points. The error 
band is even larger for the neutron data which do not participate in the fit. But within the 
la band the neutron data are well described by this fit. Due to the strong increase of the 
la band, however, we can not make precise predictions for the time-like form factors close 
to threshold. 



C. Nucleon Radii and Coupling Constants 

In Table HJ we give the nucleon radii extracted from our fits in the SC and explicit 
pQCD approaches. The first number gives the value for our best fit, while the numbers in 
parentheses indicate the range from the la band. Our values are compared to the results of 
Re f. |23l| and other recent determinations from low-momentum transfer data 5^, z2> ZH 

zUST 

The nucleon radii are generally in good agreement with other recent determinations using 
only low-momentum-transfer data given in the table. In particular, the squared neutron 
charge radius is in good agreement with the experimental value. In previous analyses |2pLE3|. 
this radius was constrained to the experimental value and not a prediction. Our result for the 
proton charge radius, however, is somewhat small. This was already the case in the earlier 
dispersion analyses of Refs. 0, HH • We speculate that the reason for this discrepancy lies 
in inconsistencies in the data sets. In this type of global analysis all four form factors are 
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3.6 3.8 4 4.2 4.4 



t [GeV ] t [GeV ] t [GeV ] 

FIG. 7: The nucleon electromagnetic form factors for time-like momentum transfer with the explicit 
pQCD continuum (right panel: G^f , left panels: Gp 1 ). The proton data participate in the fit while 
the neutron data are a genuine prediction. The solid line gives our best fit while the dashed lines 
indicate the error band obtained by the la deviation as discussed above. 





SC approach 


explicit pQCD app. 


Ref. [23] 


recent determ. 


r p E [fm] 
r p M [fm] 
(r™) 2 [fm 2 ] 
i% [fm] 


0.844 (0.840 . . . 0.852) 
0.854 (0.849 . . . 0.859) 
-0.117 (-0.11 ... -0.128) 
0.862 (0.854 . . . 0.871) 


0.830 (0.822 . . . 0.835) 
0.850 (0.843 . . . 0.852) 
-0.119 (-0.108 ... -0.13) 
0.863 (0.859 . . . 0.871) 


0.848 
0.857 
-0.12 
0.879 


0.886(15) [72, 73, 74] 
0.855(35) [Z& 25] 
-0.115(4) [52] 
0.873(11) [76] 



TABLE I: Nucleon radii extracted from our fits in the SC (2nd column) and explicit pQCD 
approaches (3rd column). The first number gives the value for our best fit, while the numbers in 
parentheses indicate the range from the la band. For comparison, we give the results of Ref. 
(4th column) and other recent determinations from low-momentum transfer data 5^, 7^, 2^> 
3,0 (5th column). 
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analyzed simultaneously and both data at small and large momentum transfers enter. This 
can be an advantage or a disadvantage depending on the question at hand. 

The discrepancy is not likely to be explained by 2j physics. In Ref. [t3] , it was shown that 
27 exchange (including intermediate nucleons only) has a very tiny effect on the extraction 
of the proton radius from ep scattering data. We have also performed various fits where the 
proton charge radius was constrained to values between 0.88 and 0.90 fm. The quality of 
these fits is not quite as good as for our best fits but still acceptable (x 2 /dof ~ 2.6). The 
fits exhibit a pronounced bump-dip structure in G P E . However, they are on the low side of 
the data in the interval Q 2 = 0.3. ..0.8 GeV 2 . These conclusions remain valid if we use a 27 
corrected data basis for the proton form factor (including intermediate nucleons in the 27 
exchange only) j78[. 

We have also extracted the ujNN couplings from our fit. The ujNN coupling constant is 
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related to the residua via 




i = 



1,2 



(24) 



where f w = 17 is the electromagnetic coupling of the u>. The u> resonance couplings, e.g., play 
an important role in addressing the issue of isospin violation in the nucleon form factors [79] . 
While the vector residua of the u are fixed relatively well by the fits, = 0.60 . . . 0.83 GeV 2 , 
even the sign of the tensor residua can not be determined, a 2 = —0.13 . . . 0.37 GeV 2 . This 
leads to the following range for the uNN coupling constants: 



We have not extracted a coupling constant for the (j) to the nucleon since the strength 
appears also in the KK and pit continua and the interpretation of the residual <fi pole 
strength is ambiguous. 

D. Inclusion of the preliminary CLAS data 

The CLAS collaboration at Jefferson Lab has recently taken new data for G n M in the 
range 0.6 < Q 2 < 5 GeV 2 [7]]]. These data are still preliminary and were not included in the 
fits discussed above. In the this subsection, we present a fit in the explicit pQCD continuum 
approach where the preliminary CLAS data are included. In Fig. |HJ we show our results for 
space-like momentum transfers compared to the published world data (blue circles) and the 
preliminary CLAS data (green triangles) [7l|. Again, we get a good description of most data 
within our error bands. In the region 0.5 GeV 2 < Q 2 < 1.0 GeV 2 the preliminary CLAS data 
differ significantly from most of the published world data. The reason for this discrepancy 
is not yet understood. Our fit prefers the CLAS data in this region. Some of the world data 
are even out of our la band. The description of the timelike data and the nucleon radii in 
this fit is of similar quality as for the fits described in subsections IV Al and IV Bl 

In addition to the lu, this fit has two more isoscalar poles (M Sl ~ 1.1 GeV and M S2 w 
1.4 GeV) and 3 isovector poles (M Vl « 1.0 GeV, M V2 w 1.6 GeV, and M„ 3 « 1.8 GeV). 
Moreover, it contains an explicit pQCD continuum term, Eq. (J2TJ) . as discussed above. 
This fit has 15 free parameters and a total x 2 /dof of 2.2. The slight increase in the total 
X 2 compared to the similar fit in subsection IV Bl is due to the inconsistency between the 
preliminary CLAS data and the older data in the region 0.5 GeV 2 < Q 2 < 1.0 GeV 2 . The 
parameters of this fit are given in detail in Tables IVI and IVT1 in the Appendix. 

E. Pion Cloud of the Nucleon 

Friedrich and Walcher (FW) recently analysed the em nucleon form factors and performed 
various phenomenological fits 27]. Their fits showed a pronounced bump-dip structure in 
Ge which they interpreted using an ansatz for the pion cloud based on the idea that the 
proton can be thought of as a virtual neutron-positively charged pion pair. They found 
a very long-range contribution to the charge distribution in the Breit frame extending out 
to about 2 fm which they attributed to the pion cloud. While naively the pion Compton 
wave length is of this size, these findings are indeed surprising if compared with the "pion 
cloud" contribution due to the 2ir continuum contribution to the isovector spectral functions 
discussed in Sec. IIII Al 



9uiNN 



16. 7.. .23.1, and 9ujnn = 



3. 6. ..10. 3 . 



(25) 
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FIG. 8: The nucleon electromagnetic form factors for space- like momentum transfer with the 
explicit pQCD continuum. The blue circles indicate the published world data while the green 
triangles give the preliminary CLAS data for G 7 ^ jzj]]. The solid line gives our best fit while the 
dashed lines indicate the error band obtained by the la deviation as discussed above. 

As was shown by Hammer, Drechsel, and Meifiner [2^], the 2tv continuum contributions 
to the long-range part of the nucleon structure are much more confined in coordinate space 
and agree well with calculations in ChP T 14711 and earlier (but less systematic) calculations 
based on chiral soliton models, see e.g. j8fj. In the dispersion-theoretical framework, the 
longest-range part of the pion cloud contribution to the nucleon form factors is given by the 
2tt continuum - the lowest-mass intermediate state including only pions. 

As a consequence, it remains to be shown how the proposed long-range pion cloud can be 
reconciled with what is known from dispersion relations and ChPT. In order to clarify this 
issue, we have performed various fits in order to understand what structures in the spectral 
function are required to reproduce the bump in G%. We find that the structure can only 
be reproduced if additional low-mass strength in the spectral function below t < 1 GeV 2 is 
allowed beyond the 2tt, KK, and pn continua and the u> pole. (See also Ref. [5^| for some 
preliminary results on this question.) Since the spectral function is well understood in this 
region in terms of meson continua and vector meson dominance, such strength was explicitly 
excluded in the fits of Sec. El 

In Fig. El we show the results in the explicit pQCD approach with additional low-mass 
strength allowed. In this fit, all constraints were removed, the neutron charge form factor 
behaviour in the region of the bump-dip structure was enhanced by artificially lowering the 
error bars on the experimental data as seen by the fit. The results for all four form factors 
in the space-like region are compared to the world data. 
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FIG. 9: The nucleon electromagnetic form factors for space-like momentum transfer with bump-dip 
structure for G\. The solid line gives our best fit while the dashed lines indicate the error band 
obtained by the la deviation as discussed above. 



In general, we get a good description of all data within our error bands. In particular, 
we obtain the desired bump-dip structure in G% at the cost of low-mass poles, which appear 
close to the u mass in the isoscalar channel, and close to the 3-pion threshold in the isovector 
channel. The latter pole is weakly coupled. The behavior of the proton charge form factor 
around Q 2 ~ 1 GeV 2 is not resolved. The fit has 22 parameters and two effective poles 
(one isoscalar and one isovector) come out below 1 GeV. In addition to the lu, this fit has 3 
isoscalar poles and 4 isovector poles. One of the isovector poles is a broad resonance with 
a width parameter of order 13 GeV. The time-like data are not included in the fit. For the 
space-like data alone the total x 2 /dof is 0.9. The parameters of this fit are listed in detail 
in Tables IV11I and IV111I in the Appendix. 

While this fit gives a good description of the space-like data, it is in contradiction to 
much of what is known about the structure of the spectral function below 1 GeV. Finally, 
we note again that the error bars of the data on G 1 ^ in the region of the bump-dip structure 
had to be lowered artificially in order to obtain the desired structure. Taking it seriously 
requires to overcome the constraints imposed by unitarity and analyticity in the spectral 
function below 1 GeV. 
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VI. SUMMARY & OUTLOOK 



Dispersion theory simultaneously describes all four nucleon form factors over the whole 
range of momentum transfers in both the space-like and time-like regions. It allows for the 
inclusion of constraints from other physical processes, unitarity, and ChPT and therefore is 
an ideal tool to analyze the form factor data. 

We have presented the results of our new form factor analysis. The spectral function has 
been improved compared to earlier analysis in various respects, (i) It contains the updated 
2n continuum [38j], as well the KK 3|| 36[ and pn continua j3?J as independent input, (ii) 
The pQCD behavior of the form factors at large momentum transfer has been included in 
two different ways: using super convergence relations and a broad resonance to mimick the 
QCD continuum (SC) and by including an explicit pQCD continuum term. Moreover, we 
have generated la error bands for all our fits by performing a Monte Carlo sampling of all 
fits with a total x 2 /dof in the interval [xlnim Xmin + 1-04]. 

Our fits give a consistent description of the world data in the space-like and time-like 
regions. We find that some residual <fi pole strength is still required in addition to the 
KK and p7r continua. The FENICE data for the neutron time-like form factors do not 
participate in the fit but are reproduced within the errors bands. This is an improvement 
over our earlier studies 21, 2^]. The nucleon radii are generally in good agreement with 



other determinations with the exception of the proton charge radius which comes out smaller 
for our best fits. As discussed in subsection IV CI enforcing a larger radius r p E ~ 0.88 fm 
leeds to acceptable fits with a slightly higher \ 2 . While the ujNN vector coupling constant 
is determined relatively well by the fits, even the sign of the tensor coupling constant can 
not be determined. We have also performed fits including the preliminary data for G n M from 
the CLAS collaboration at Jefferson Lab [71]. For lower Q 2 these data are in conflict with 
earlier data, still our fits seem to prefer the trend set by the CLAS data. The bump-dip 
structure in G\ advocated by Friedrich and Walcher 27( can only be obtained by artificially 
lowering the error bars of the data and allowing additional strength in the spectral function 
below 1 GeV. 

There are two important extensions of this work. First, one could fit directly to the 
Coulomb-corrected cross section data, which would eliminate possible inconsistencies in the 
form factor data basis [75|. This could be further refined by consistently removing the full 27 
exchange from the cross section data. The latter, however, remains a formidable challenge 
because of the contribution of the nucleon excited states. 
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APPENDIX A: FIT PARAMETERS 



In this Appendix, we give the values of the fit parameters for the various fits. Table |H] 
contains the parameter values for the fit in the SC approach. 



Resonance 


TV It T rl 

Mass [GeV] 


a\ [GeV J 


a 2 [GeV z J 


r [GeV] 


U 


0.782 


0.755960 


0.370592 




<f> 


1.019 


-0.776537 


-2.913229 




Sl 


1.124860 


0.902379 


2.484859 




S2 


2.019536 


0.022798 


-0.130622 


5.158635 


Vl 


1.062128 


-0.127290 


-2.162533 




V 2 


1.300946 


-1.243412 


3.704233 




v 3 


1.493630 


4.191380 


-7.091021 






1.668522 


-3.176013 


3.723858 




V5 


2.915451 


0.048987 


0.075965 


19.088297 



TABLE II: Fit parameters in the SC approach. This fit has 17 free parameters and a total x 2 /dof 
of 1.8. 

In Table IIHI we list the resonance parameters for the fit in the explicit pQCD approach. 
The parameters for the explicit pQCD term are given in Table ITVl 



Resonance 


Mass [GeV] 


oi [GeV 2 ] 


a 2 [GeV 2 ] 


UJ 


0.782 


0.616384 


0.114681 





1.019 


0.159562 


-0.329255 


Sl 


1.799639 


0.128654 


0.026174 


Vl 


1.000000 


-0.309199 


-1.078695 


V2 


1.627379 


3.695960 


-4.301057 


V3 


1.779245 


-3.693109 


3.630255 



TABLE III: Resonance parameters for the fit with explicit pQCD continuum. This fit has 14 free 
parameters and a total x 2 /dof of 2.0. 





a\ 


h [GeV" 2 ] 


ci [GeV" 1 ] 


0.002321 


-0.028391 


0.152903 


0.161871 


a s 2 


al 


b 2 [GeV" 3 ] 


c 2 [GeV" 1 ] 


-0.126598 


-0.011693 


1.159998 


1.150000 



TABLE IV: Parameters of the explicit pQCD term for the fit with explicit pQCD continuum. 
This fit has 14 free parameters and a total x 2 /dof of 2.0. 

In Table E] we list the resonance parameters for the fit including the preliminary data for 
G n M from the CLAS collaboration at Jefferson Lab Yl] in the range 0.6 GeV 2 < Q 2 < 5 GeV 2 . 
The parameters for the explicit pQCD term for the fit including the preliminary CLAS data 
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Resonance 


Mass [GeV] 

L J 


ay [GeV 2 ] 


a 2 [GeV 2 ] 

L J 




0.782 


669166 


—0 135957 




1 045277 


-0 025807 


0.001144 


S2 


1.400423 


0.261240 


-0.053588 


Vl 


1.022008 


-0.279441 


-1.215307 


V2 


1.644552 


3.823047 


-4.561225 


V3 


1.770845 


-3.849954 


4.027035 



TABLE V: Resonance parameters for the fit including the preliminary CLAS data for 
This fit has 15 free parameters and a total x 2 /dof of 2.2. 



71] 



a\ 


a\ 


by [GeV" 2 ] 


c x [GeV" 1 ] 


-0.000186 


-0.026941 


0.219241 


0.169695 


a s 2 


a\ 


b 2 [GeV" 3 ] 


c 2 [GeV" 1 ] 


0.000527 


-0.001835 


0.004155 


0.106343 



TABLE VI: Parameters of the explicit pQCD term for the fit including the preliminary CLAS 
data |7jj for G 1 ^. This fit has 15 free parameters and a total x 2 /dof of 2.2. 



for G 1 ^ are given in Table IVTl 

In Table IVH| we list the resonance parameters for the fit with the bump-dip structure 
and additional low-mass strength allowed. The parameters for the explicit pQCD term in 



Resonance 


Mass [GeV] 


ay [GeV 2 ] 


a 2 [GeV 2 ] 


T [GeV] 


u 


0.782 


-3.088199 


1.516336 




sy 


1.087524 


-9.309347 


5.152311 




S2 


0.857075 


7.969599 


-3.102716 




S3 


1.145783 


5.332548 


-3.754331 




Vl 


0.315028 


0.002785 


-0.008642 




V2 


1.523890 


-3.257202 


3.767630 




v 3 


1.323997 


2.770486 


-5.497436 




V<y 


2.834388 


0.177584 


-0.011050 


13.477161 



TABLE VII: Resonance parameters for the fit with the bump-dip structure. This fit has 22 free 
parameters and a total x 2 /dof of 0.9 (space-like data only). Note that S2 and v\ are the additional 
low-mass poles necessary to generate the bump-dip structure. 



this fit are given in Table IVIII1 
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af 


a\ 


b x [GeV" 2 ] 


c a [GeV" 1 ] 


-0.786259 


-0.320320 


0.971368 


1.235451 


a| 


al 


b 2 [GeV" 3 ] 


c 2 [GeV" 1 ] 


-0.000484 


0.033410 


0.091209 


0.994702 



TABLE VIII: Parameters of the explicit pQCD term for the fit with bump-dip structure. This fit 
has 22 free parameters and a total x 2 /dof of 0.9 (space-like data only). 
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